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ABSTRACT 

Recent work by Levitan et al has expanded the long-term photometric database for AM CVn stars. In par¬ 
ticular, their outburst properties are well-correlated with orbital period, and allow constraints to be placed 
on the secular mass transfer rate between secondary and primary if one adopts the disk instability model 
for the outbursts. We use the observed range of outbursting behavior for AM CVn systems as a func¬ 
tion of orbital period to place a constraint on mass transfer rate versus orbital period. We infer a rate 
x 1O _9 M 0 yr -1 (P or b/1000 s) -5 2 . We show the functional form so obtained is consistent with the 
recurrence time-orbital period relation found by Levitan et al using a simple theory for the recurrence time. 
Also, we predict their steep dependence of outburst duration on orbital period will flatten considerably once 
the longer orbital period systems have more complete observations. 

Subject headings: accretion, accretion disks - binaries: close - novae, cataclysmic variables - stars: individual 
(AM Canum Venaticorum) 

1. INTRODUCTION 

Cataclysmic variables (CVs) are semi-detached binaries 
consisting of interacting stars in which a Roche-lobe filling 
secondary transfers matter to a more massive, and also more 
compact primary. The dwarf novae (DNe) constitute a sub¬ 
class further characterized by semiregular outbursts of several 
magnitudes, recurring on timescale of days to years (Warner 
1995a). The novalikes reside at higher mass transfer rates and 
do not show DN outbursts. The limit cycle accretion disk in¬ 
stability model (DIM) has been successful in explaining the 
outbursts (Smak 1984). CVs evolve to shorter orbital periods, 
driven by angular momentum loss (AML) from a combination 
of a magnetic wind from the secondary star and gravitational 
radiation (Knigge, Baraffe, & Patterson 2011 = KBP). Al¬ 
though the DIM is generally agreed to be the correct explana¬ 
tion for DN outbursts, mainly because the observed dividing 
line between steady and outbursting systems (i.e., DNe and 
novalikes) agrees with theory (Smak 1983b), further progress 
has come slowly. For instance, any theory for the recurrence 
time for DN outbursts (Cannizzo, Shafter, Wheeler 1988 = 

CSW) is muddied by the fact that the outburst properties at a 
given orbital period exhibit a wide scatter (Warner 1995a). 

The AM CVn stars are a subset of CVs at very short or¬ 
bital period, less than an hour, whose spectra are dominated 
by helium (Nelemans 2005, Solheim 2010). The prototype, 

AM CVn, was discovered neary fifty years ago (Smak 1967). 

Paczyriski (1967) proposed the system to be a short orbital 
period binary with two degenerate, He-rich stars with a pe¬ 
riod evolution driven by gravitational wave radiation (GWR), 
but he considered detached rather than semi-detached bina¬ 
ries. Motivated by observational inferences of mass transfer 
in AM CVn, namely rapid photometric flickering (Warner & 
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Robinson 1972), Faulkner et al. (1972) presented the first 
self-consistent model for AM CVns, correctly taking into ac¬ 
count their semi-detached nature. 

The coincidence in stable periodic variations in spectral 
lines with photometric variations can sometimes be used to 
obtain orbital periods, and the weak X-ray emission indicates 
a white dwarf (WD) primary accretor rather than a neutron 
star (Nather et al. 1981, Patterson et al. 1993, Groot et al. 
2001). AM CVn stars are thought to have evolved beyond the 
minimum orbital period which divides nondegenerate and de¬ 
generate secondaries. The systems span a wide range in mass 
transfer rate Mt . In fact, in the AM CVn stars one sees a 
range spanning not only the high to intermediate range equiv¬ 
alent to the novalike —>■ DNe transition in normal hydrogen- 
containing CVs, but also a range spanning intermediate to low 
Mt values. For the low Mt systems the disk is too cool to 
have DIM outbursts. 

In some sense the AM CVn stars are a better laboratory for 
the DIM than normal DNe because less scatter is expected in 
their Mt values at a given orbital period, given that GWR 
is the sole AML mechanism. KBP find that the normal (i.e., 
solar composition) short orbital period DNe below the period 
gap must have AML enhanced by a factor ~2.47 above that 
given solely by GWR in order to produce the observed period 
gap. This only applies to an ensemble average; as with the 
DNe above the gap, the outburst properties of systems below 
the gap show considerable scatter at a given orbital period. In 
addition, Mt varies by more than four orders of magnitude 
across the range of orbital periods in the AM CVn systems, 
which makes orbital period-dependent properties of the out¬ 
bursts more noticeable. 

Amassing a large database for AM CVn stars has been ham¬ 
pered by their faintness. Now, thanks to the results of dedi¬ 
cated transient surveys, an avalanche of new systems and data 
on outbursting behavior has become available. Levitan et al. 
(2015) present the results of a comprehensive study of AM 
CVn systems over nearly 10 yr. They present outburst data 
on 32 systems with known orbital periods, ranging from 5 to 
65 min. A similar study with a slightly smaller sample was 
presented by Ramsay et al. (2012). 
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In Section 2 we discuss the DIM in the context of AM CVn 
accretion disks. In Section 3 we look at the recurrence time 
for superoutbursts, and in section 4 the superoutburst dura¬ 
tion. Section 5 summarizes our findings. 


2. THE DIM IN AM CVN SYSTEMS 
2.1. Background 

A calculation of the vertical structure of geometrically thin 
disks reveals the steady state physics underlying the limit cy¬ 
cle behavior, the vertical structure as effective temperature 
T e g versus surface density E (Meyer & Meyer-Hofmeister 
1981, Cannizzo & Wheeler 1984). One finds a hysteretic rela¬ 
tion between midplane temperature T m ;d and surface density 
E, or equivalently, between effective temperature T e g and E. 
Each radius r and viscosity parameter a value has its own S- 
curve. A parameterization of the results of these calculations 
allows a determination of the surface densities at the max¬ 
ima and minima in the S-curve, E max (r, a) and E m ; n (r, a), 
as well as other physical quantities associated with these ex¬ 
trema, like midplane and effective temperature. One relates 
the local accretion rate to T e g jmax and T e g jln i n using the stan¬ 
dard Shakura & Sunyaev (1973) relation 


crT eS A 
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( 1 ) 


In the DIM, gas accumulates in quiescence and accretes 
onto the central object in outburst, (e.g., Cannizzo 1993a, 
Lasota 2001 for reviews). The phases of quiescence and out¬ 
burst are mediated by the action of heating and cooling fronts 
that transverse the disk and bring about phase transitions be¬ 
tween low and high states, consisting of neutral and ionized 
gas, respectively. During quiescence, when the surface den¬ 
sity E(r) at some radius within the disk exceeds a critical 
value E max (r), a transition to the high state is initiated; dur¬ 
ing outburst, when E(r) drops below a different critical value 
S m i n (r), a transition to the low state is initiated. Low-t>high 
transitions can begin at any radius, whereas high—dow transi¬ 
tions begin at the outer disk edge. 

This situation comes about because of the following: In 
quiescence the disk is very non-steady so that mass accumu¬ 
lates in the outer regions. The surface density distribution 
is bounded by E max (r), which increases with radius. (Both 
Emax(t) and E m i n (r) scale close to linearly with r.) In the 
outburst disk, however, E(r) oc r -3 / 4 . The accretion disk 
mass is conserved, i.e., the disk mass accumulated by the 
end of quiescence AM co id is the same as that in the hot disk 
A A/hot immediately after the heating transition has occurred. 
Therefore there must be a substantial redistribution of E(r) 
- from a profile oc r in quiescence to oc r -3 / 4 in outburst. 
Since E(r) is smallest at large radii in the outbursting disk, 
and since E(r) < E m i n (r) is the condition for the cooling 
transition to begin, cooling fronts are always initiated in the 
outer disk. 

One can define a “maximum mass” AA/ nlax = 
f 2nrdrT, max (r) that the disk could possibly achieve during 
the accumulation phase, i.e., quiescence. Obviously one can¬ 
not have E(r) > E max (r) in quiescence or else the instability 
would have already been triggered. In practice, time depen¬ 
dent calculations show that the true disk mass at the end of 
quiescence AM co id is typically ~1/10 — 1/3 of the “maxi¬ 
mum mass” (e.g., Cannizzo 1993b). Therefore one may write 
AM cold = /AM max , with / ~ 1/10 - 1/3. 


Levitan et al. (2015) restrict their attention to superout¬ 
bursts in their outburst statistics. These are analogous to 
long outbursts in DNe above the period gap (van Paradijs 
1983); most of the mass accumulated during quiescence ac¬ 
creted onto the primary before the disk shuts off. For nor¬ 
mal, “short” outbursts, only a few percent of the stored gas 
accretes onto the central object: the thermal time scale of thin 
disk is short compared to the viscous time scale, and the cool¬ 
ing front that is launched from the outer edge of the disk al¬ 
most as soon as the disk enters into outburst traverses the disk 
and reverts it back to quiescence. For disks that have been 
“filled" to a higher level with respect to AM max , the surface 
density in the outer disk can significantly exceed the critical 
surface density E m ; n . In order for the cooling front to be¬ 
gin, however, the outer surface density E(r ou ter) must drop 
below E m i n (r ou ter)- Disks in this state generate much longer 
outbursts, with slower “viscous” plateaus, because the entire 
disk must remain in its high, completely ionized state until 
enough mass has been lost onto the primary for the condition 
S(r ou ter) < £ min (r 0 uter) t0 be satisfied. 

Various studies have investigated the DIM in AM CVn sys¬ 
tems (Smak 1983a, Cannizzo 1984, Tsugawa & Osaki 1997 = 
T097, El-Khoury, & Wickramasinghe 2000, Menou, Perna, 
& Hernquist 2002, Lasota, Dubus, & Kruk 2008 = LDK, 
Kotko et al. 2012 = KLDH). These investigations generally 
consist of first calculating the steady state accretion disk struc¬ 
ture by integrating the vertical structure equations to param¬ 
eterize the S-curve relation between surface density E(r, a) 
and effective temperature T e ff(r, a), and then using scalings 
for the steady state physics as input into a time dependent 
model to calculate light curves. T097 did not solve the full 
set of equations for the vertical structure but rather prescribed 
a functional form for the flux F(z). We restrict our consider¬ 
ation of the scalings to the two most recent studies LDK and 
KLDH, which integrate the complete set of structure equa¬ 
tions and present a complete set of scalings for both the local 
minima and maxima in E. 

As a model for the mass-losing secondaries in AM CVns, 
Deloye et al. (2005) calculate pseudo-evolutionary sequences 
for donors with varying degrees of degeneracy. In their Fig¬ 
ure 2 they calculate tracks in the Mx-Porb plane for four 
isotherms based on central donor temperatures T c ranging 
from 10 4 K to 10' K, assuming a constant primary mass 
0.6 Mq. They also indicate the upper and lower bounds of 
the instability strip for the DIM taken from T097. For their 
two lowest tracks, T c = 10 4 K and 10 6 K, the instability strip 
spans roughly the correct (i.e., observed) period range. Their 
tracks do not represent true evolutionary sequences since T c 
is taken to be constant along a track. 

2.2. Instability Criteria as Power Law Scalings 

The range of mass transfer from the secondary star feeding 
into the outer disk Mt which allows for unstable behavior, 
i.e., dwarf nova outbursts, is set by the local stability criteria 
at the inner and outer edge of the accretion disk. 

KLDH calculate many S-curves for accretion disks relevant 
for ultracompact binaries with no hydrogen, X = 0. The ba¬ 
sic finding is that the steady state scalings for the DIM are 
shifted to higher surface densities and temperatures. They 
present three sets of scalings for (i) Y = 1, (ii) Y = 0.98, 
Z = 0.02, and (iii) Y = 0.96, Z = 0.04. Their results 
are given in terms of a, r, and M\. The range for unsta¬ 
ble disk behavior is determined by the S-curve for the in- 
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ner and outer disk radii. If Mj > Mr,2 = M m - m , outer, 
the rate of accretion associated with the minimum in £ at 
fouter, the disk will be stable in the high, ionized state. If 
Mt < Mt,i = Afmax, inner- the rate of accretion associated 
with the maximum in £ at ri nner , the disk will be stable in the 
low state. 

The KLDH scalings are not convenient; for comparison 
with observations we must replace disk radius r with orbital 
period P or b in the KLDH scaling for Mt, 2 , and with the pri¬ 
mary radius in the Mt, 1 scaling. We follow T097 in making 
these conversions: 

( 1 ) For the outer scaling, Mt, 2, we relate the outer disk ra¬ 
dius r outer to orbital period using Figure 5 from van Haaften 
et al. (2012). We identify the Roche lobe radius Rli with 
r outer and fit a simple power-law to Rli/ a as a function 
of the mass ratio q = m 2 /mi 0 The semi-major axis is a. 
Over the relevant range 0.01 < q < 0.03 we fit Rli/cl ~ 
0.85(<7/0.01) -0 06 . Based on the results of van Haaften et al. 
( 2012 ) we adopt a fiducial secondary star 77 i 2 (P or b) relation 

"’ 2 = 0038 G®ft) 3 <2) 

(see their Sect. 2.3), corresponding to R 2 05 M/ -0 ’ 18 , which 
gives, using Kepler’s law. 



where r o , 10 = r oute r/10 10 cm. 

(2) For the inner scaling, Mt, 1 , we adopt the standard inner 
boundary zero torque condition (Shakura & Sunyaev 1973) 
for which the maximum effective temperature 

= 0.488 < 4 ' 

is reached at 49/36 times the inner edge, taken to be the pri¬ 
mary radius. (The “max” in this equation refers to the local 
radial maximum in T e g- in a steady state disk due to the in¬ 
ner boundary condition, not the local maximum in £ in the 
S-curve.) Thus Mt,i is enhanced by a factor (0.488) —4 and 
evaluated at (49/36)P w d. 

For the primary mass-radius relation, we fit a power law 
to Eggleton’s scaling of a zero-temperature WD, i.e., for 
which the ratio of atomic number to atomic weight Z/A = 
1/2 (Rappaport et al. 1987, see their eqn. [19], with 
Mch = 1.44 Mq). Taking the scaling relevant for a He 
WD (van Haaften et al. 2012, see their eqn. [25], with 
M p = 5.66 x 10” 4 M 0 ), we tabulate values for x = log JV%d 
and y = log 7 ?wd over the range of interest, 0.5 < m\ < 0.7, 
and fit a least squares power law 7 ?wd = 10 880 cm mi -0 ' 62 . 
Over the fit range this relation gives a maximum deviation 
< 1.5% from the Eggleton scaling. 

Applying these conversions to the scalings given in KLDH 
(and ignoring the weak a dependencies) we obtain scalings 
for Mt i and Mt 2 relevant for their three compositions (i) 
Y = 1 ,’ (ii) Y = 6.98, Z = 0.02, and (iii) Y = 0.96, Z = 

5 We use Mi and M 2 to refer to the primary and secondary masses in cgs 
units, and mi and m 2 to indicate M/Mq. 


TABLE 1 

KLDH-based coefficients for local extrema 


comp. 

Mt, 1,0 
(M 0 yr -1 ) 

Mt, 2,0 
(M 0 yr“i) 

e lml 

£2ml 


e P 

(i) 

10-10.85 

10 —8-60 

- 2.54 

0.16 

0.89 

2.0 

(ii) 

to - it- 04 

10-8.82 

-2.50 

0.16 

0.89 

1.99 

(iii) 

lo-n 10 

10 -8.93 

-2.49 

0.16 

0.88 

1.97 


0.04. Adopting the general forms 

Mt, 1 = M t , 1,0 TOi(F2) £lml ( 5 ) 

Mt ,2 = M t , 2,0 rm(Pi) e2ml (1 + q ) e * ( ' , ( 6 ) 

\1000 s J 

we may now calculate the coefficients. These are given in 
Table 1. 

Levitan et al. (2015) presents a list of AM CVn systems 
ordered by P or b, with outburst properties indicated. The rate 
of mass transfer Mt decreases sharply with P or b, so that the 
shortest P or b systems have disks in permanent outburst, those 
with intermediate P or b exhibit outbursts, and those with the 
longest P or b have disks in permanent low states. The divid¬ 
ing point between high state systems and outbursting systems 
lies at Pi « 20 min. The dividing point between outbursting 
systems and low state systems, P 2 , is not as straightforward. 
A block of systems starting at P or b =44.3 min are listed as 
not showing outbursts 0 , but two systems among these do show 
outbursts, at P or b = 47.3 and 48.3 min. The expressions for 
Mt, 1 indicate a steep inverse scaling with mi, therefore it 
seems probable that these systems have somewhat high pri¬ 
mary masses. Therefore we set P 2 = 44 min. Furthermore 
we adopt roi(Pi) = mi(P 2 ) = 0.6. We note that the nil 
values for AM CVn systems are not well-constrained obser- 
vationally from dynamical measurements. 

We may now use the observed instability strip for AM CVn 
systems to constrain the secondary mass transfer rate Mt- 
Let us assume a power law Mt = A (P or b/1000 s) n . For 
simplicity we adopt the convention that Mt > 0. Since we 
expect Mt to decrease with orbital period, we set Mt,i = 
Mt{P 2 ) and Mt ,2 = Mt(Pi), which gives two equations in 
two unknowns. We may solve for the normalization constant 
A and dependence on orbital period n = din Mt /dinP or b- 
This yields the general solution 

log[{*T,2,o/*T,i.o) 9 (A/1000 s)‘>] 

n = - SEaTa)- ,7) 

4 = M Mmi (A)-- (AMY’, (» 

where 

g = (l + g) e «m^i- £ i m i. (9) 

Table 2 gives the values for A and n for the three KLDH 
compositions, adopting Pi = 20 min, P 2 = 44 min, and 
mi = 0.6 ± 0.05. The magnitude of the putative assigned 
error on mi was propagated through to n and A in order to 
indicate the strength of the dependency. The main uncertainty 
entering into n and A is the assumed primary mass. For in- 

6 There is also a system within the instability zone at 35.2 min indicated as 
not having outbursts. We disregard it on the assumption it must be anomalous 
in some way, or its outbursts may have been missed. 
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3. RECURRENCE TIME FOR OUTBURSTS 


TABLE 2 

Coefficients for M t = A(P orb /1000 s) n with mi = 0.6 ± 0.05 


composition 

A(M 0 yr- d ) 

n 

(i) 

8.82 ±0.6 x 10" 9 

-5.38 ±0.3 

(ii) 

5.27 ±0.4 x 10 -9 

-5.23 ±0.3 

(iii) 

3.88 ±0.3 x 10" 9 

-5.06 ±0.3 


stance, taking mi = 1 would give n — —7, a much steeper 
ATr(-Porb) relation. 

Since precise abundances for AM CVns as a group are 
not known, we adopt a representative fiducial scaling Mr — 
5 x 1O _9 M 0 yr _1 (P orb /1000 s) -5 ' 2 , relevant for the middle 
range of scalings given in KLDH. 


2.3. Considerations from AM CVn Binary Evolution 

How does this Mt law compare to theoretical expectations? 
As mentioned above, in AM CVn systems the driving of the 
mass transfer is most likely the angular momentum loss due 
to GWR. For that case one can derive the expected scaling 
of mass transfer rate with orbital period as is done in Warner 
(1995b). We begin by expressing the mass transfer rate in 
terms of the angular momentum loss 


AM 7 ) oc MMM 1 + M 2 I m 

M i (-Crb/owR a 

(e.g., see Savonije et al. 1986; Marsh et al. 2004; eqn. [6] of 
van Haaften et al 2012), where M 2 = —Mt- Assuming Mi 
remains constant and M2 -C Mi and using Kepler’s law to 
replace orbital separation with orbital period, we find 


M2 

M 2 


oc M 2 P-1 /3 . 


( 11 ) 


We then can use the fact that the size of the Roche lobe for 
given period depends only very weakly on the mass of the ac- 
cretor (the well known period-mean density relation) to find 
how P orb scales with M2, given the mass-radius relation for 
the donor R 2 oc M|. Hence 

Porb OC (X M 2 (3C “ 1)/2 (12) 

and therefore 

M 2 oc M|p-® /3 oc P 0 4 / b (3C_1) “ 8/3 - (13) 

For ( = —1/3 we find an exponent n = —14/3. In van 
Haaften et al. (2012) the values of ( are plotted in their Fig¬ 
ure 2. For low masses ( > —1/3, yielding a larger absolute 
value of the exponent. Indeed, the fit they make to the depen¬ 
dence of mass transfer rate on orbital period (their Appendix 
A) gives n = —5.32 for a 1.4 Mq accretor. Warner (1995b) 
uses ( = —0.19, based on Savonije et al. (1986), and finds 
n = —5.21. We conclude that the exponent for the scaling we 
infer from the DIM based on the observed instability strip for 
AM CVns is in good agreement with expectations from stellar 
structure if we adopt mi ~ 0.6. 


7 Note that there is an exponent “—1” missing in eqn. (4) of Warner 
(1995b); for £ = —1/3 the donor mass and orbital period scale inversely, 
not linearly. 


For 11 of the 32 AM CVn systems with known orbital peri¬ 
ods in their study, Levitan et al (2015) have enough coverage 
to make quantitative statements about their outburst proper¬ 
ties. They find a relation for the recurrence time for outbursts 

/ p \ 7 - 35 

frecur = 1-46 d —— + 24.7 d. (14) 

recur ^ 1000 sj 

They only include superoutbursts. Their relation (see their 
Figure 12a) is much tighter than the comparable plot for DNe 
above the minimum period (e.g., see Figures 2-4 of CSW; 
Figure 11 of Patterson 2011). 

What is the expectation from DIM for the recurrence time 
and is it consistent with this scaling? 

CSW formulated an analytical expression for t le cur . Their 
full complexity is not needed, and indeed one can take a rather 
simple approach. Regardless of whether one subscribes to the 
thermal-tidal instability for superoutbursts or the plain DIM 
(Osaki & Kato 2013), the normal outbursts in a system ex¬ 
hibiting both normal outbursts and superoutbursts represent a 
tiny fraction of the mass budget. Therefore a good approxi¬ 
mation is that superoutbursts are the only outbursts, and that 
during a superoutburst essentially all the mass stored in the 
cold state is accreted. The recurrence time is then simply 

_ AM C Qid _ / AM max _ / / 2'itrdrYi max (j "5 

t'recur = ~ r ~ = = Tw ■ 

Mt Mt Mt 

(15) 

where AM max is the maximum mass that could be stored in 
quiescence and / is the fraction the disk is filled, relative to 
this maximum. 

For specificity we adopt the middle of the three compo¬ 
sitions considered by KLDH. Their £ max for Y = 0.98, 
Z = 0.02 is 

£ max = 612 g cm’ 2 a c _r°- 82 no 1 ' 10 mf 0 ' 37 , (16) 

where a c -i is the value of a in the low state of the disk, a co id, 
normalized to 0.1. Hence 

AM max = 1.04 x lO" 9 M 0 ctc-i -0 ' 82 mi 0 85 

(1+9) 1 - 03 (17) 

v \ 1.000 s ) ’ 

Evaluating the recurrence time gives 

frecur = 7.59 d /_! « c _i' 0 ' 82 TOi°- 85 (1 + g) 103 

(18) 

= 4 - 92d (it®/) 751 (19) 

for mi = 0.6, where /_i = //0.1 = a c _i = 1. Thus for 
a mass transfer rate with n = din Mt/ d\n P or b = —5.2, 
the simplest recurrence time scaling gives t lecur oc / :> orh 7 ' 5 , 
which is close to the observed relation. 


4. DURATION TIME FOR OUTBURSTS 


For the AM CVn superoutburst durations Levitan et al. 
(2015) find 


fdur = 0.89 d 


p 


4.54 


10.6 d. 


1000 s 


( 20 ) 
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This relation appears to be less reliable than their f recur (-Porb) 
relation: for the four longest period systems only upper limits 
are given. In fact, for the other six systems, those with 22 min 
< Porb < 29 min, their data are consistent with Tiur(-Porb) 
being constant with orbital period. We note that the four out¬ 
burst durations given as upper limits in Table 2 of Levitan et 
al. are plotted as actual values in their Figure 12, panel 3, 
and enter into their power law fitting. As noted in their Ap¬ 
pendix A, for the four systems with only upper limits (UL); 
on fdur i they arbitrarily take fdur i = 0.75(UL)i, with error 
0.25 (UL)i. 

What would one expect from theory? 

There are several ways to construct an outburst duration 
timescale for superoutbursts. During a long outburst with a 
“viscous plateau”, the cooling front cannot propagate due to 
excess surface density at the outer disk edge, relative to E m ; n , 
therefore the only option is for accretion onto the primary so 
that S(r ou ter) is gradually reduced. Therefore it is reasonable 
to consider the viscous time scale in the outer disk in the hot 
state as an approximation to the superoutburst duration. 

We follow the method given in King & Pringle (2009). 
They estimate a local peak accretion rate during outburst as 
~ 2AM hot /fvi S c, where AAfhot is the mass of the hot disk. 
The outbursting disk mass is determined by that stored in qui¬ 
escence, AAfhot = AAfcoid- This eliminates M from the 
standard equations given in Frank, King, & Raine (2002), 
which relate the locally defined viscous time to a, r, and Af. 
Using the outer disk edge for r and relating it to orbital period 
as previously gives 

f v i sc = 15.1da^ 1 8 rn;- 34 (l + 5 ) 0 ' 16 < 21 > 

where a h -i = a h ot/ 0 . 1 . 

This is much flatter than the Levitan et al. relation for 
Liur(-P or b), but their fitted relation is dominated by includ¬ 
ing upper limits for (for systems with P or b > 30 min) as 
part of their fit, which makes a direct comparison with theory 
problematic. 

5. DISCUSSION AND CONCLUSION 

Previous workers have examined the stability properties of 
AM CVn systems in the context of expectations from the 
DIM, e.g., T097 (see their Fig. 4), Nelemans (2005, see his 
Figs. 1 and 6 ), LDK (see their Fig. 3), and KLDH (see their 
Fig. 3). The last work, KLDH, has provided the most com¬ 
plete study to date. Although they do not provide explicit 
formulae between the critical mass transfer rates and orbital 
period as in T097, KLDH plot them in their Fig. 3, as do the 
aforementioned studies. Our main difference with KLDH is 
that for Mt,i we also take into account the zero torque bound¬ 
ary condition at the inner edge of the disk (following T097), 
in addition to just considering the primary mass. However, 
Mt,i depends steeply on the primary mass (oc mj~ 2 ' 5 ), and 
this uncertainty will likely dominate that associated with the 
boundary condition refinement. KLDH plot in their Fig. 3 
the expected secular Mt versus P or b from stellar evolution¬ 
ary models, and find consistency with the DIM in terms of 
stability of observed systems. Thus the difference between 
our study and KLDH is one of perspective: KLDH combine 
observations with evolutionary models to show that AM CVn 
outbursts can be explained by the DIM, whereas we assume 
that the DIM is correct in order to derive Mr (Porb ), and then 


explore the ramifications of the derived law vis a vis not only 
A/rf/'orb) predicted from evolutionary models, but also the 
resultant outburst properties versus P or b- 

In summary, we apply the DIM model to the recent re¬ 
sults on AM CVn systems obtained by Levitan et al. (2015), 
using scalings for the helium-rich accretion disks in such 
systems taken from Kotko et al. (2012). The orbital peri¬ 
ods defining the edges of the instability strip Pi and P 2 per¬ 
mit us to infer a mean secondary mass transfer rate Mr — 
5 x 10~ 9 Mq yr _1 (P or b/1000 s) -5 2 . Our finding of a steep 
inverse dependence n ~ —5.2 is consistent with theoretical 
expectations, but our result is dependent on taking mi ~ 0.6; 
higher m\ steepens the relation. Treating mi as a variable 
and all other parameters on the right hand side of equation (7) 
as constant we may write n ~ —5.23 — 7.75 log(mi/0.6) for 
Y = 0.98, Z = 0.02. The largest uncertainty in n enters via 
mi. 

We emphasize that the power-law form for P or b is not an 
outcome of our analysis but an assumption, valid in our ap¬ 
plication only over 20 min < P or b < 45 min. However, 
our inferred P or b values at 20 and 45 min are in line with 
those estimated for stable systems with similar orbital periods 
(KLDH, see their Fig. 3). We note that the precise values 
of P 0 rb defining the edge of the instability strip. Pi and P 2 , 
also affect n and A, although not as strongly as mi. Lastly, n 
and A are only weakly dependent on the composition of the 
gas and the viscosity parameter a; these two factors enter via 
A7 t.2,o/4Tt.i,o (see eqn. [7]). This ratio is relatively insensi¬ 
tive to composition and virtually independent of a. The full 
dependencies of n and A on all input parameters are given in 
eqs. (7)-(9). 

The simplest possible theoretical expression for the re¬ 
currence time from the DIM gives t recur oc r^ ter / Mr oc 
Porb 2 ' 31 -^^ 1 , so that (dlnt recur /<ilnP or b) = 2.31 — 
(—5.2) ~ 7.5, close to the value 7.35 found by Levitan et 
al. (2015). However, a larger assumed value of mi in our 
model would increase the exponent. Thus, if the DIM is rele¬ 
vant for AM CVn outbursts, the primaries must have masses 
mi ~ 0.6. The Levitan et al. constraint on the outburst dura¬ 
tion oc P 4 r b 4 appears to be dominated by incompleteness for 
the upper half of their range in orbital period. Our theoretical 
prediction is that a more complete time sampling of AM CVn 
outbursts, especially at longer orbital period, will ultimately 
reveal a much flatter law fdur 9S Poib- 

We thank Joe Patterson for organizing a small CV work¬ 
shop at Columbia University in Fall of 2014 which provided 
a stimulating environment for discussion. We also thank the 
anonymous referee whose comments improved the paper. 
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